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In this article the effect of noise on Grover’s algorithm is analyzed, modeled as a total depolarizing 
channel (TDCh) and a local depolarizing channel in each qubit (LDCh). The focus was not in error 
correction (e.g. by the fault-tolerant method), but to provide an insight to the kind of error, or 
degradation, that needs to be corrected. In the last years analytical results regarding mainly the 
TDCh model have been obtained. In this paper we extend these previous results to the local case, 
concluding that the degradation of Grover’s algorithm with the latter is worse than the former. It 
has been shown that for both cases with an TV-dependent small enough error-width, smaller than 
1/VN for total error and l/(-\/TVlog 2 N) for the local case, correction is not needed. 
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I. INTRODUCTION 


Quantum systems cannot be completely isolated from the environment. When a quantum system is to be controlled 
externally, e.g. gate applications, state preparations and others, additional errors are generated. Such effect on 
quantum algorithms has been studied lately by several authors |l|-|5|. 

While there is a fault-tolerant model of quantum computing based on the correction of errors below a certain 
threshold 00 , this method is very expensive in computational resources. Therefore, it is important to diagnose the 
performance of quantum algorithms affected by noise before making any correction. 

Grover’s quantum algorithm is an oracle-based search of an unordered database that, without noise, presents a 
quadratic speedup over the brute-force classical search algorithm |M3. However, this advantage is affected drastically 
if the oracle is faulty 01 or in the presence of noise El- 

Local error models are important due to the fact that any implementation of a quantum circuit are affected by gate 
errors. In this article we study the effect of a known rate of noise on Grover’s search algorithm, extending the work 
done in El with total depolarizing channel error model (TDCh), and also analyzing the effects of local depolarizing 
channel errors (LDCh). For both model errors, it is shown that quantum speedup disappears when the error rate 
is constant. Furthermore, we have concluded that the degradation of Grover’s algorithm with local error (LDCh) is 
worse than with total error (TDCh). This is coherent with recent results EH- 

The article is divided as follows: in order to introduce the equations used later, in section 2 Grover’s algorithm 
is briefly explained. Sections 3 and 4 deal with the effect of total and local depolarizing error (TDCh and LDCh) 
on Grover’s search, respectively. Some conclusions are drawn in section 5. The major calculus are given in the 
appendixes. 


II. GROVER’S QUANTUM SEARCH ALGORITHM 


Grover’s quantum search algorithm is known to be optimal, in some sense, to solve the problem of finding a marked 
element in a unsorted database of N elements using kc r = f \/ivj oracle queries 0 . Such problem may be shortly 
stated as follows: suppose we have a database of N = 2 n of quantum states, being n the number of qubits that span 
the Hilbert space ( Ti. = C 2 ), and an unknown marked state among them. Given an oracle ( or black box ) to identify 
whether an element is the one being searched, the goal is to find such marked state with high probability and in as 
few steps as possible. 

Let \t) denote the target basis state and |s) = - 7 = X^o 1 1*) the superposition of all basis states. To be able to 

introduce quantum error, we use the density matrix notation (H = C 2 x C 2 ), obtaining: pt = |t)(f| and po = |s)(s|, 
respectively. The algorithm is represented in Algorithm |T| 


Algorithm 1: Grover’s search algorithm 

1. Set up the superposition state po- 

2 . Apply the oracle operator O = 2\t) (t| — I. 

3. Apply the diffusion operator D = 2 |s) (s| — I. 

4. Repeat steps 2 and 3 ^ n/aJ — 1 times. 

5. Perform measurements in the canonical basis in each qubit. The target state will emerge with high probability as 
1V» 1. 


It is straightforward that after applying the oracle operator and diffusion operator k times, the result is 


p(k) = G k P o(G^) k , (1) 

where G = DO is also known as the Grover operator. It can be shown that the density operator obtained afterwards 
is 


p{k) = |sfc)(s fe |, 


( 2 ) 


where 
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|sfc) = sin ((2k + 1 ) 0 ) \t) + cos(( 2 k + l) 0 )|f) 

N -1 


l*> = 


1 


i^t 


|f) and 8 = 


arcsm 


1 


Hence, the probability of success (measuring the marked element) after k steps is 

p{k) = sin 2 (( 2 k + 1 ) 0 ). 


(3) 


(4) 


III. GROVER’S SEARCH ALGORITHM WITH TOTAL DEPOLARIZING CHANNEL 


In this section we analyze how the error modeled as the Total Depolarizing Channel (TDCh) affects Grover’s search 
algorithm. Similar work, but with different focus, has been done by Vrana et al Fl~4| . Instead of minimizing the mean 
cost as done in fbij we propose a maximization of the probability, because a closed form expression can be obtained. 
This enables us to compare the results with the error modeled as the Local Depolarizing Channel (LDCh), in Sec. 
ED The TDCh is an error model, which maps an n-qubit state p to the maximally mixed state jt with probability 7 , 
and leaves it unchanged with probability (1 — 7 ) 17[, given by 


e(p, 7) = (1 - l)p + l~^- 


(5) 


This can be interpreted in the probabilistic sense as a Binomial distribution, where success or failure should be 
interpreted as having p or respectively. We will use these terms indistinctly. 


A. State evolution with error 

We now consider the effects of the TDCh in Grover’s algorithm. That is, every step of the algorithm involves 
applying the Grover operator, followed by the TDCh error. Since the TDCh error commutes with any unitary 
operator, the density matrix obtained after k steps is 

P(k, 7) = (1 ^ 7 ) k p(k) + (!-(!- 7 )*)^, (6) 

where p(k) is given by Eq. ([3]). This has a very straightforward interpretation: we have probability of success (1 — y) fc 
and probability of failure 1 — (1 — 7 )^. 

This implies the probability of finding the marked element is 

m 7 ) = (1 - 7 ) k p(k) + 1 ~ ( 1 A T 7)fc , (7) 

where p(k) is given by Eq. (jTj. Notice how the maximally mixed state jt introduces uniform probability over all 
possible states. 

As can be seen from Fig. [I] and Fig. [51 the first maxima of the probability p(k), at step fc max , moves left with 
respect to kc r as either 7 or N increase [3|. This can also be appreciated on Fig. [3] which also shows that for any 
number of qubits, e.g. n = 10, the behavior is similar. A good approximation of k max can be deduced by maximizing 
Eq. ® (see®, resulting in 


fc ma x( 7 ) = m ax 


7 r — arcsin <5 — arcsin ([l — -A] S) 
40 


( 8 ) 


where 


6 = 


1 



4 9 

111(1-7) 


2 


(9) 
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FIG. 1. Probability vs. number of steps with n = 10 qubits and different values of 7 : 0 (A), l/(4%/]V) (B), 1/%/A (C), 4/%/iV 
(D) and 1 (E). 
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FIG. 2. Probability vs. normalized number of steps with 7 = 0.01 and different number of qubits n: 10 (A), 12 (B), 14 (C) 
and 16 (D). 


and 9 is given by Eq. fl3J. We provide two approximations of Eq. ©: for n 1 (N = 2 n ) and 0 < 7 <C 4= 


^max (ff) 


7 ry/N V 7 

— 8~~ 


and for n !§> 1 with 1 > 7 » —7= 

7TV N 


^max (7) 


max 


1 

In (1 - 7 ) 


1 

2 




( 10 ) 


( 11 ) 


Notice how for Eq. (USD. the approximation is valid only for small values of 7 ; while the approximation given by 
Eq. (1111) is valid for a large range of 7 . Fig. [I] shows different approximations of fc max . 
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FIG. 3. Number of steps fc max /fcG r vs. width of noise 7 vs. number of qubits n. 
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FIG. 4. Step of maximum probability vs. width of noise 7 , with n = 10 qubits and different approximations of fc max : obtained 
by simulation (A) and as in equations (0 (B), (flOl) (C) and (fill) (D). 


B. Order estimation of the algorithm 


If we apply classical amplification and stop at step fc, we have the order of the algorithm given by the estimate 
mean cost (expected number of oracle queries): 

MC Kl := MC[k,p(k, 7 )] = —(12) 

P(k, 7 ) 

It is well-known that the classical brute-force search algorithm has order Q(N) with a constant factor of 1/2. 


1. Stopping the algorithm at ltGr 

If we were to stop the algorithm at we would have the following estimate cost for the algorithm: 

_ kcr 

7 (1 - 7 ) kGr p (k Gr ) + j? (1 - (1 - r y) kGr )' 


MC kG . 


(13) 
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Since we consider JV > 1, we have: p (£y?r) ~ 1 and kc r = 

MC kGrtl ~- 


: J 


VN ~ ttVn 


, leaving Eq. m as: 


t y/~N 


(1-7 ) S2 P + jf (l - (1-7)^) 

Notice that for any fixed 7 > 0 the order is O (Ny/N). For small values of 7 ( 2L ^7 <C 1) Eq. (THU) becomes 


(14) 


, ^ „ 7 xy/N ( 7 t\/N 

M Ck Gr)7 ~ — "7 — 1 4 Z - 7 


(15) 


It is interesting to find the width of noise that keeps the quantum speedup (order of the algorithm Q(\^N)). Hence, 
it is straightforward that 7 needs to be a function of N. In fact, one needs 7 -C 1/kcr, as shown in 14]. 

However, Grover’s quantum search algorithm can still outperform the brute-force classical search algorithm. In 
fact, there is a maximum width of noise (qdassicai) that leaves the mean cost of Grover’s quantum search algorithm 
Y- Such is found by using Eq. 


/classical 


/—1 H- 2 /cgv \ kGr 

V -l + N ) 


which for N 1, becomes 


'"/classical 


Un(AS) 

7 Ty/N 


(16) 


(17) 


2. Stopping the algorithm at k ma , x 


If we were to stop the algorithm at fc ma x( 7 ), Eq. (fl2l) becomes MC [fc m ax(7),p(^max(7))7)]- Therefore, in contrast 
to stopping the algorithm at kc r , for any fixed 7 > 0 the algorithm is of order <d(N). For large values of N (N 1) 
and small values of 7 (v/iVy/ (2tt) <C 1 ), the mean cost becomes 


MC k 


k(7)>7 


tVN 


tVn 


1 - 



(18) 


Comparing this result with Eq. m , we can see that an extra constant factor of approximately 0.80 appears. Then, 
we need 7 <C 1 /( 0 . 8 &Gr) to mantain the quadratic speedup (order of the algorithm Q(y/N)). 

For any fixed 7 > 0, the mean cost becomes 


MCu 


1 ( 7).7 


N 

9 — 87 


(19) 


As a result, the maximum width of noise that leaves MC [fc max ( 7 ),p(^max( 7 ), 7 )] classical (N/2) is independent of N, 
be. 7classical ~ 7/8. 


IV. GROVER’S SEARCH ALGORITHM WITH LOCAL DEPOLARIZING CHANNEL 

The Local Depolarizing Channel (LDCh), maps an n-qubit state p to a mixed state applying the depolarizing 
channel in every qubit independently, i.e. 

e(p) = £i(p, a) o £ 2 (p, a) o ■ ■ ■ o e „(p, a), (20) 

where £j(p, a) corresponds to the depolarizing channel (Eq. (f3|l) acting in the qubit i, and a the probability of error. 
Analogously to Sec. IIII1 we define the probability of obtaining the marked element after k steps as p L (k,a), and the 
step at which the first maxima of the probability occurs as /c^ ax (o:). 
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A. Operator-sum representation for Grover state evolution with LDCh 

For one qubit states the Depolarizing Channel error (Eq. (0) can be expressed using Kraus’ operators (Operator- 
sum representation), as 

4 

s(p,a)=Y, M ipM}, ( 21 ) 

1=0 


where M 0 = yj (l — x) a ° an d = (* e {1, 2, 3}), being a M the Pauli matrices. 

In the case applying local errors in all qubits in a n-qubits state there are 4" operators (all posible combinations), 
which are of the form 0 



where cr,; 9 is the i q Pauli matrix (i q £ 0,..., 3) and m is the amount of er 0 (/) operators present in a particular M c . 

If we consider one step of the algorithm as applying the Grover operator G and the LDCh, we get the following 
state after k steps 


4 n 4 n 

p L (k, a) = [° n Pini °n 

Cl Ck 


(23) 


where 


k 

O u = l[(M cl G ), 
1=1 


(24) 


l represents the step, each M Cl has the form of Eq. (E^l) and p m i = |s)(s| is the initial Grover state. 
For the first step (k = 1), such probability is 


which for n» 1 becomes 


p L (l,a) 


1 


2 n (2 ri 


23n—4 

+ — (2 n - 4 f 
16 v ; 



p L (l,a) 



(25) 


(26) 


Due to the exponentially increasing complexity of Eq. (l23l) . we propose to analyze p L (k,a) using lower and 
upper bounds in terms of the TDCh model (because of simplicity and relatively tight bounds), and a first order 
approximation. 


B. Probability Bounds 


Grover’s algorithm with LDCh error can be interpreted as a ternary tree where each depth-level corresponds to a 
step. The initial node is given by pi n i = |s)(s|, and in each node we have the paths shown in Fig. [5] In fact, the 
actual probability of getting the marked element at a step k is obtained by taking all the possible paths (which are 
4 nk if one uses the Kraus’ operators). 

In order to find upper and lower probability bounds, we assume the best and worst case for the desired probability, 
by taking ^unknown = GpG ^ or p U nknown = //AT, respectively. Thus, the probability bounds are 


Pu(k,a) 


Pi(k, a) 


(1 - a n ) k p(k) + 
(1 - a) nk p(k) + 


1 - (1 - a n ) k 


N 

1 (1 a ) nk 


(27) 


N 


(28) 
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Punknown 

FIG. 5. Ternary tree at step k of the algorithm with LDCh. 


that hold 


Pi(k, a) < p L {k, a) < p u (k, a), (29) 

Vfc € N, a £ [0,1]. Since these are of the form of the TDCh (Eq. ([ 6 ])), we define these bounds as their equivalent 
width of error 7 (a), i.e. 

7«( a ) = <*”, (30) 

= 1 — (1 — a) n . (31) 

Even though 7 i(a) provides a reasonable probability bound, 7 u (a) does not (see Fig. O- A better approach is to 
replace it by 


1u{ol) 


na 

2 + na 


(32) 


C. First order approximation of the LDCh 

Here we will propose a method for analyzing the effect of the LDCh in Grover’s algorithm in terms of its first order 
approximation (in a, as in 0 ), proving also that this approximation serves both as an approximation for small values 
of width of error a and as a lower probability bound for any a. 

After applying the LDCh in every qubit and every step (Eq. (l23l) h the probability of success after k steps is of the 
form 


nk / o \ nk—i 

(*>“)(i~ t) (f)/<(«>*) 

2=0 ^ ' 


(33) 


where i is the total number of errors up to step k. Hence, the first order approximation is 


p L (k, a) 


1 - 


3a 


nk 


fo(n,k) + 1 - 


3a 


nk —1 


(|) h(n,k), 


(34) 


where fo{n, k) = p(k) is Grover’s algorithm original probability and /i(n, k ) corresponds to the all the contributions 
of applying one error in any qubit at any step of the algorithm (see[B|). As can be seen in Fig. [Gl Eq. (13411 is also a 
lower bound of the probability because all fi(n,k ) are non-negative (sums of probabilities). 
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(a)Probability evaluated at kc r vs. width of error a for n = 8 qubits and its bounds. 
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(b)Detail for small values of a. 


FIG. 6 . Comparison between different bounds of probabilities: numerical simulation (A), first order approximation (B), lower 
bound with 7 i(ot) = 1 — (1 — a) n (C), upper bound with 7 u (a) = na /(2 + net) (D) and upper bound with 7 u (a) = a n (E). 


D. Order estimation of the algorithm 

Here we denote the estimate mean cost of the algorithm with LDCh error by MCfc a := MC L [fc,p L (/c, a)]. Using 
the monotonicity of both p L (k, a) and ( 7 / — y u ) in a, and Eq. (l29l) , we have found bounds for the mean cost as follows: 

MC k , 7u <MC£ a <MC k!7l . (35) 

That is, given a fixed number of steps k, the bounds p(k,j u ) and p(fc, 7 ;) induce bounds on the mean cost. 


1. Stopping the algorithm at kGr 


Considering kc r as the final step, we have an estimate mean cost of the algorithm given by MCfc G a , bounded as 
in Eq. (l35l) . Similar to the TDCli case, for any fixed a > 0 the order is Q(Ny/N). 
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For N 1 and small values of a ( akc r log 2 IVC 1) the mean cost becomes 

MC^ a = ^ + (aNn 2 log 2 N , (36) 

where £ is a constant in the interval [1/32,1/16] (from bounds). Using the first order approximation, Eq. (l34l) and 
/i(n, fc), we have found that ( is approximately 70/2048. The latter is much closer to 1/32 than to 1/16, and hence, 
the actual probability is much closer when using y u than when using ji. Because of the simplicity and the validity of 
the expression, we encourage to use j u to estimate the probability p L (k, a). 

From Eq. (l36l) . one can see that to mantain the quadratic speedup a <C 2/{kGr log 2 N) is needed. 


2. Stopping the algorithm at fc max 


We now consider the induced bounds for /c max , which are justified by the continuity of p(fc, 7 ) in both 7 and n. 
These are 


l max. I — a max — A, max_u 5 


(37) 


where fc max _; = fcmax ( 7 z) and fc max .„ = fc max ( 7 „). Since we are stopping at the step of maximum probability (in each 
curve), the following bound holds: 


MCu 


< MC^_ a < MC kn 


(38) 


In fact, using Eqs. (l38l) and (l26l) . one finds that k^ ax —► 1 and the mean cost is given by l/p L (l, a). Therefore, again, 
for any fixed a > 0 the order is Q(N). 

For N 1 and small values of a (0.55aJVlog 2 N <C 1), the mean cost is approximately 


MC^ a = ^ + C aN log 2 N, (39) 

where C is a constant in the interval [(-7T 2 — 4) /32, (7r 2 — l) /16]. And for the algorithm to be of order 0(\//V), one 
needs a 2.5 /(kcr log 2 N). 

This is coherent with the results presented with the TDCh model, where for small values of 7 , choosing kc r or & max 
does not change the order of the algorithm, but for any fixed 7 > 0 , it changes drastically. 

Analogously to the work done in Sec. IIIII one can find classical: the maximum value of a that leaves the mean cost 
N/2. Considering the worst case yields 


^classical 


< 


14 

log 2 N ' 


In contrast to the TDCh case, this width of error does depend on N. 


(40) 


V. CONCLUSIONS 

In this article we have studied the effect of noise in Grover’s quantum search algorithm, based in two models of error: 
total depolarizing channel (TDCh) and local depolarizing channel (LDCh). Our focus was not in error correction (e.g. 
by the fault-tolerant method), but to provide an insight to the kind of error, or degradation, that needs to be corrected. 
Our work extends the results found by Vrana et. al. regarding only the TDCh model, concluding that the LDCh 
error model is more devastating than the TDCh error. We also show that, for both cases, with an V-dependent small 
enough error-width, correction is not needed. 

The devastating effect of the LDCh was, in some sense, expected: in the case of TDCh model the state is partially 
altered, but the rest remains intact (the state evolves partially in the plane formed by |s) and 1 1), or otherwise becomes 
I/N). Whereas, with an LDCh error model, the state is mostly removed from the plane formed by |s) and 1 1). 

We have found a closed form for the step at which the maximum of probability occurs for the TDCh error model 
(fcmax), and provided bounds for the LDCh error model (based on the TDCh model). This step could be used as an 
indicator of the impact of the error in the algorithm. This means, comparing fc max with fcg r one has an idea of how 
degraded the algorithm is. Analogous comparisons could be done with the probability at such steps, but are, in some 
sense, equivalent. 
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From this study we conclude that to maintain the original speed up of Grover’s algorithm, different requirements 
are needed depending on the error model: in the TDCh case, the error-width has to be smaller than 1 /y/N (as in 
0 ); and for the LDCh model, it needs to be smaller than 1 / (y/N \og 2 N). 

Furthermore, there are widths of error for both the TDCh and the LDCh model for which the algorithm still 
outperforms the classical brute-force algorithm. It is interesting to note that for the TDCh model this width does not 
depend on N , i.e. for widths below 7/8 the algorithm is better than the classical (Eq. (fTT)l) ). Nevertheless, for the 
LDCh model, this classical threshold error a c i ass i ca i, depends of N (Eq. (l40l) l and goes to zero with V, reinforcing 
the idea that the LDCh degrades more the algorithm than the TDCh. 

The depolarizing channel is a particular error model. It would be interesting to study other error models taking into 
account experimental implementations of the algorithm. Future work should be done with this focus, and possibly 
including error correction. 


Appendix A: Finding fc n 


We wish to find fc max which maximizes Eq. 0 - We will assume k £ R, and then perform some adjustments. 
Considering p(k) as in Eq. 0, we find the first derivative of Eq. 0 


dp(k, 7 ) 
dk 


= (! - 7 ) 


29 sin (2(2 k + 1)9) + ( p(k) — — ) In (1 — 7 ) 


(Al) 


Assuming 7 / 1 and nulling the derivative yields 


40sin ((2 k + 1)0) cos ((2 k + 1)0) + |^sin 2 ((2 k + 1)0) — — j In (1 — 7) = 0. (A2) 

The former can be expressed as 

— C = — A sin 2 a; + H sin a; cos a:, (A3) 

where x = (2k + 1)0, A = — In (1 — 7 ), B = 40 and C = — In (1 — 7 )/Ab This can be reduced to 

— z = sin (x) cos (x + 4>), (A4) 

where z = C/y/A 2 + B 2 and (j) = arcsin (A/y/ A 2 + B 2 ). Considering x > 0 Eq. sm has two family of solutions 

X\ = 7T m — — arcsin(2^: — sin ((/)) — and 

X 2 = — (2ttiti + arcsin( 2 z — sin(</)) — (j) + n) , (A5) 

being m an arbitrary integer. Thus, the solution to the original problem is 

tt — 29 — arcsin ( 6 ) — arcsin ([l — b) 


40 


where 


5 = 


\ 1+ 

Adjusting border values with 7 —> 0 and 7 —> 1, and taking integer part yields 

7 T — arcsin S — arcsin ([l — <5) 


(A 6 ) 


(A7) 


fcmax ( 7 ) = m & x 


40 


(A 8 ) 
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Appendix B: Finding fi(n,k) 


As stated before, fi(n, k) corresponds to the depolarizing channel acting in any qubit at any step of the algorithm. 
Because of the simmetry of Grover’s algorithm in the qubits, one can assume without loss of generality that the last 
qubit is the affected one and then multiply by the number of possible ways of having the depolarizing channel acting 
on one qubit ( n ). We can also assume, without loss of generality, that the target state is \t) = |00 • • ■ 0) = |0). 

The depolarizing channel acting on one qubit (without the identity operator and considering that error has been 
commited) can be expressed with the operators 

X := /®" _1 ® <r x , Y := I®"- 1 ® a y , Z := 7®"" 1 ® <r„ (Bl) 

where / represents the 1-qubit identity operator and a M the Pauli matrices. 

It is useful to see that the set {|0), |s), |1), |p}} is closed under the operations of {X, Y, Z, G}, where G is Grover’s 
operator and | p) is defined as the normalized superposition of the even states of the canonical base as 


\p) 



(B2) 


T 

The state obtained after k steps is represented as <f>(k) = [a* bk Ck dk} , which corresponds to the quantum 
state | tpk) = «fe|0) +6fc|s) +c*,|l) + dk\p}- The probability of obtaining the marked state (|0)) is given by [tr (G.</>(fc))] 2 3 4 5 6 , 
where C = [l A/2 0 A/\/2] and A = 2/y/N. Hence, the desired function is 

k 3 

fi(n,k) = EE [tr(CA k ~ l B j A l </>( 0))] 2 , (B3) 

l=i j =l 


where (j>{k) = [010 0] T is the initial step of the algorithm (| s)) and A, Bi,B 2 , B 3 correspond to applying G, X , Y, 
or Z respectively (up to a global phase). Such matrices are given by 


A = 


-1 -A 0 -v^A 
A A 2 — 1 —A i/2(A 2 - 1) 


Hi = 


0 0 10 
010/2 
10 0 0 


0 

0 

0 

1 


0 0 

0 

-1 


'0 

0 

-1 0 


' 1 

0 

0 

o' 

b 2 = 

0 

1 

0 V2 

and B 3 = 

0 

-1 

0 

0 

1 

0 

0 0 

0 

0 

-1 

0 


0 

-y/2 

0 -1 


0 

-y/2 

0 

1 


(B4) 


(B5) 


The expression of /i(n, k), derived here, has been used to perform the series expansion made in Eq. (l36l) . The full 
expression is too large to include in this article. 
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